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Abstract — This work proposes the application of fractal descriptors to the 
i Oh analysis of nanoscale materials under different experimental conditions. We 

obtain descriptors for images from the sample applying a multiscale trans- 
form to the calculation of fractal dimension of a surface map of such image. 
Particularly, we have used the Bouligand-Minkowski fractal dimension. We 
applied these descriptors to discriminate between two titanium oxide films 
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prepared under different experimental conditions. Results demonstrate the 
discrimination power of proposed descriptors in such kind of application. 

{Sj Keywords — Porous Titanium Oxide, Nanomaterials, Pattern Recogni- 

tion, Fractal Dimension, Fractal Descriptors. 
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1 Introduction 

Morphological material characterization is a very important and challeng- 
ing task throughout Material Science, since such properties are determinant 
in the suitability of a given material to a specific application. Among sev- 
eral materials studied recently, a great deal attention has been given to 
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photoactive ones, such as titanium oxide. Titanium oxide films prepared 
electrochemically can show a wide variation in their morphology depending 
on the experimental conditions in which they are prepared. As an example, 
the morphology can range from self-organized nanoporous (9 19 21 up to 
nanostructures without pores definition [3}[5]. As described, many authors 
have shown that morphology is an important aspect to be considered in 
the photoactive properties of Ti02 films (3j[l3] . Recently we have shown 
that the photoactivity of Ti02 films prepared by galvanostatic anodization 
is strongly affected by the morphology on the early stages of the anodization 



process 18 . In [18] we have investigated the effect of morphology using an 
Image Analysis Method, which is based on the computational analysis of 
visual attributes, such as color, shape and texture. Among these attributes, 
texture is a powerful characterizer for material analysis. Although the con- 
cept of texture has no precise definition, it may be comprehended as the 
spatial organization of pixels in a digital image. A physical consequence of 
this definition is that this attribute is capable of express characteristics such 
as luminosity and roughness of a digitalized object. In this way, it allows a 
robust description and discrimination of material images, particularly Ti02 
samples used in this study, which show few or no information about the pore 
diameter, but have a rich morphology related to texture information. 

Unlike conventional textural data, natural textures do not present any 
evident quasi-periodic structure, but persistent random patterns [11] . Among 
the approaches employed for texture analysis, fractal methods are proposed 
as the best solution [12| . These methods, based on fractal dimension [14] 



or multifractal spectrum 10 , can measure the complexity of an object tex- 
ture, which corresponds to the shape irregularity, also related to the spatial 
occupation of an object. Therefore, fractal measurements are capable of 
quantifying the texture homogeneity, allowing a comparison among their in- 
formation and the consequent discrimination of original materials. Using as 



example the materials here investigated, some authors [16 , 20 have shown 



that the fractal dimension of TiC>2 films can also be correlated to photoac- 
tivity properties of these materials. However, in that papers, the analysis is 
valid for very regular morphologies, showing deviations when more complex 
surfaces are considered. 

Although the fractal dimension is a good descriptor to characterize a 
texture image, it is inefficient in applications involving the discrimination of 
a large amount of objects. In fact, it is easy to find objects with different 
aspects presenting the same fractal dimension [12| . In order to solve this 
drawback, the literature presents approaches which extract a lot of descrip- 
tors based on fractal geometry, such as Multifractals [To] and Multiscale 
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Fractal Dimension (MFD) [4j[8j[T5]. In the literature, Backes et al. [T] ap- 
plied MFD in texture analysis employing volumetric Bouligand-Minkowski 
MFD (VBMFD) to extract a set of descriptors from natural textures with 
very good results. VBMFD is based on the intensity image mapping onto 
a 3D surface. In the following, this surface is dilated by a variable radius r 
and the volume V{r) is calculated for each radius. This expansion process 
gives a precise measure of the pixel arrangement. As the radius r grows, 
an interaction among dilation spheres is observed, interfering in V(r) value. 
Therefore, these values capture changes in spatial distribution of textures 
along different scales. In this way, VBMFD uses V(r) values as descriptors 
for a texture, capable of distinguishing different textures with their different 
spatial arrangement. 

Considering the exposed above, in this work it is proposed the use of 
VBMFD descriptors to discriminate among nanoscale Ti02 films prepared 
electrochemically using two different experimental conditions and charac- 
terized by field emission gun scanning electron microscopy (FEGSEM) . 
In order to demonstrate the efficiency of the proposed method, this work 
obtains VBMFD descriptors from these films. 

2 Materials and Methods 
2.1 Fractal Theory 

Since its introduction by Mandelbrot [14], the literature presents a lot of 
works using fractal theory to describe and discriminate several kinds of ma- 
terials [16||20| . Most of them use fractal dimension as a descriptor of the 
original samples. This is explained by the fact that fractal dimension mea- 
sures the complexity of the object, related to the irregularity or the spatial 
occupation of that object. This property is strongly correlated to physical 
important properties. 

2.1.1 Fractal Dimension 

The fractal dimension Y\M is the most commonly used measure to charac- 
terize a fractal object. Despite its importance, we cannot find a unique 
definition for this concept. The most ancient definition corresponds to the 
Hausdoff-Besicovitch dimension. 

If X € is a geometrical set of points, the Hausdoff-Besicovitch di- 
mension dim,H(X) is given by: 

dim H (X) = inf {s : H S {X) = 0} = sup {H S {X) = oo} , (1) 
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where H S {X) is the s-dimensional Hausdorff measure, given by: 

H S (X) = lim inf j \Ui\ s : Ui is a 5-cover of X | , (2) 



. i=l 



where || expresses a diameter in K n , that is, |J7| = sup\x — y\ : x,y G U. 

Otherwise, the Hausdorff-Besicovitch definition application could be im- 
practicable in many real situations. This is the case of discrete objects 
represented in a digital image, as those studied here. For such applications, 
we have an alternative definition of fractal dimension, which is a general- 
ization from the topological dimension. Thus, the fractal dimension D is 
provided by: 

DW-totStWi, (3) 
e->0 log(±) 

where iV(e) is the number of objects with linear size e required to cover the 
whole object X j7]. More generically, N(e) may be considered a measure 
which varies according to the scale e. Such measure is characterized by a 
power-law relation with the scale [17) . A lot of fractal dimension methods 
were developed by using distinct measures, such as Bouligand-Minkowski 
[7], box-counting [7), Fourier [17]. Here, we are focused on the Bouligand- 
Minkowski fractal dimension. 

2.1.2 Bouligand-Minkowski 

The original definition of Bouligand-Minkowski fractal dimension dim^(X) 
depends on a symmetrical structuring element Y: 

dim M (X, Y) = inf {r, meas M (X, Y, r) = 0} , (4) 

where measM is the Bouligand-Minkowski measure: 

, YV , v V(dXerY) 

meas M {X, Y, r) = hm — , (5) 

r— >o r n T 

where r is the radius of Y and V is the volume of the dilation between Y 
and the boundary dX of X. 

By using neighborhood techniques we may find a simplified version of 
this dimension, which eliminates the explicit dependence from Y: 

dim M (X) = lim (dt- WP® 1 ® ). (6) 
e^o \ loge J 



4 



For the cases such as those studied here, when X 6 5ft 3 , the topological 
dimension is Dt = 3 and Y e is a sphere with diameter e. For each value of e, 
each point in X is dilated by Y e and the number of points inside the dilated 
structure is the dilation volume V. 



2.2 Proposed Method 

A significative drawback of fractal dimension is that it is only a unique real 
number. Therefore, it is possible to find a great number of objects which 
despite its same fractal dimension present appearance completely diverse. 
This is exemplified in the Figure [l] using textures from Brodatz data set (2J. 



The authors in 15 proposed the Multiscale Fractal Dimension technique 
to solve this issue. In MFD, the fractal dimension is calculated for the 
object observed under different spatial scales and each value is used as a 
descriptor for the object. In this work, we propose the analysis of nanoscale 
FEG images by applying a multiscale approach to the Bouligand-Minkowski 
fractal dimension. 




Initially, we map the image Img £ [1 : M] X [1 : N] — > M onto a 
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threedimensional surface 

Surf = {i,j,f(i,j)\(i,j) G [1 : M] x [1 : N]}, (7) 

such that: 

fihj) = {M, -,max.gray}\f = Img(i,j), (8) 

where max-gray is the maximum pixel intensity. 

In the following, the surface is dilated by a variable radius r, as depicted 
in Figure [4j Then the dilation volume V(r) for each dilation radius is cal- 
culated. The value of V(r) also corresponds to the number of points with a 
distance at most r from the object. Therefore, the exact Euclidean distance 
transform (EDT) [6] becomes a powerful and efficient tool for this calculus. 

In the threedimensional space, EDT may be defined as the distance of 
each point in the space to a subset of it. In our case, this subset is the 
surface and the EDT for each point outside Surf is given by: 

EDT(p) = min{d(p, q)\q £ Surf}, (9) 

where d is the Euclidean distance. 

In particular, in exact EDT, the distances present discrete values E: 

E = 0,l,V2,...,l,..., (10) 

where 

i eD = {d \ d = { e+j 2 ) l i 2 -i,jen} (ii) 

The dilation volume is given through: 

r 

V{r) = Y,Qii), (12) 

i=i 

where 

Q(r) = (x,y,z)\g k (P) - [g r (P) n U^g^P)], (13) 

such that: 

g r (P) = {(x, y, z)\[(x - P x ) 2 + (y- P y ) 2 + (z - P z ) 2 } , (14) 

where 

P = (x,y,z)\f(x,y,z) £ Surf (15) 

The authors in |lj propose the use of values V(r) as descriptors for tex- 
ture. Such technique is named Volumetric Bouligand-Minkowski Multiscale 
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Fractal Dimension (VBMFD). In that work, the authors apply this technique 
to the analysis of plant leaves, achieving interesting results. It is important 
to stress out that V(r) is directly related to the Bouligand-Minkowski di- 
mension for maximum radius r. Besides, each radius corresponds to an 
observation scale, from further (greater radius) to closer (smaller radius). 

This work proposes the application of VBMFD descriptors for the anal- 
ysis of two Ti02 samples prepared under different experimental conditions. 
The goal is to discriminate between the samples based on images extracted 
from different regions of each film. 

2.3 Sample Preparation 

The samples were prepared as follows. Titanium sheet (Alfa Aesar, 99.99%, 
0.25mm thick) with an exposed area of 1 cm 2 , and two platinum sheets were 
used as working and counter electrodes, respectively. Before the anodization 
process, the working electrode was polished with #1000 SiC and then with 
#1200 SiC emery paper followed by vigorous washing with deionized water. 

The experiments were carried out under galvanostatic conditions with 
a home-made current source, measuring the potential difference between 
working electrode and counter electrode with an HP 34410A multimeter 
coupled to a computer by an in-house software routine developed with HP- 
VEE 5.0 software. The two samples investigated were prepared in oxalic acid 
(sample 1 in a concentration of 0.05 mol L , and sample 2 in 0.5 mol L^ 1 ) 
applying a constant current density of 10 mA cm~ 2 and 20 mA cm~ 2 (sample 
1 and sample 2, respectively) on the working electrode. The temperature 
was kept constant at 10 °C during the anodization of both samples. 

After the preparation, the samples were morphlogically characterized 
using a field emission gun scanning electron microscope (Supra35 - Zeiss). 
From both samples we have extracted ten images from different regions of 
the material in order to collect a representative amount of data. 

3 Results 

The fractal descriptors were used for the discrimination of two Ti02 films 
prepared under different experimental conditions. From each sample, we 
have extracted two micrographs, represented in the Figure [2] 

Initially, we apply the transform which maps the intensity image onto a 
surface, as illustrated in the Figure [3j 

Next, the surface is dilated by a variable radius r, as depicted in the 
Figure |4j The conventional Bouligand-Minkowski dimension technique is 
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Figure 2: Texture images extracted from titanium oxide under two different 
conditions (above and below). 




Figure 3: Texture mapped onto a surface, (a) Original texture, (b) Three- 
dimensional surface. 

extracted from the curve log(V(r)) x r, in which, V(r) is the diation volume. 
Here, the descriptors obtained by the VBMFD technique are provided by: 

D'(r) = log(V(r)) (16) 

In order to eliminate redundancies from descriptors D' k (r) of the k th 
image the following operation was performed: 

n 

D k (r)=D' k {r)-Y J D' l {r)/n, (17) 
i=l 

where n is the total number of images in the dataset and D k (r) is the used 
descriptor. 

Figure [5] shows the descriptor curves for each image from TiC>2 samples. 
Each line aspect (solid and dashed) corresponds to a sample prepared in a 
different experimental condition. 
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Figure 4: Dilated surfaces with different radii, (a) Original surface, (b) 
Radius 2. (c) Radius 5. (d) Radius 10. 

It is noticeable that descriptors curve aspects discriminate strongly the 
materials. We do not observe any interlacement among curves from different 
conditions. The whole graphic reflects a high inter-class variability and a 
low intra-class variation in descriptors. Such result is in accordance to the 
pixel arrangement nature, measured through fractal descriptors. 

Moreover, we notice the complex shape of each curve, reflecting the 
richness of information expressed in the descriptors. We observe a clear 
improve in the description compared to the classical fractal dimension. 

Figure [6] shows the distribution of two statistical important scores of frac- 
tal descriptors: from Principal Component Analysis (PCA) and Canonical 
Correlation Analysis (CCA). 

While PCA plot measures the dispersion of whole dataset, CCA mea- 
sures the correlation among data from a same material. Both graphics show 
a cluster of points (corresponding to the 2 main scores) in each material con- 
ditions, comprobing again the discrimination power of fractal descriptors. 
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Figure 5: Descriptors curves from each texture image. Solid curves corre- 
spond to sample 1. Dashed ones correspond to sample 2. 
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Figure 6: Distribution of correlation scores, (a) Principal Component Anal- 
ysis, (b) Canonical Correlation Analysis. 



4 Conclusions 

In this work, we proposed the application of fractal descriptors to the dis- 
crimination of materials under different experimental conditions. We ob- 
tained the descriptors by applying a multiscale approach to the fractal di- 
mension estimation of texture images extracted from the material. 
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The curves of descriptors showed a high accuracy in the discrimination. 
This result demonstrated the power of the proposed technique in the mod- 
elling of materials. It also comprobed the validity of materials analysis 
through computational texture analysis. 
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